SeCU>^tTY  CLASSlPlCATjON  OP  This  PAGC  0*i«  Enftmd) 


REPORT  DOCUMENTATION  PAGE 


\  RCPOAT  number 


i2.  GOVT  ACCESSION  NO.i 


NRL  Memorandum  Report  5155 


4.  title  ('Mtf  Su6r/r/«j 


PLASMA  EFFECTS  ON  LEVEL  SHIFTS  OF 
H  AND  He"^ 


7.  AuTHOPf*; 


M.  Blaha*  and  J.  Davis 


t.  PER^OMmING  organization  name  AnO  aooress 

Naval  Research  Laboratory 
Washington,  DC  20375 


n.  CONTROLLING  OFFICE  NAME  ANO  AOORCSS 


Office  of  Naval  Research 
Arlington,  VA  22217 


14.  monitoring  AGCkCy  name  4  aOORESV//  dtiUtfit  from  Conif^tllnt  Oitte0) 


i«.  Distribution  statement  (at  thu  Rapoet) 


Approved  for  public  release;  distribution  unlimited. 


17.  distribution  statement  (ai  IA«  in  Black  JO,  If  dificrcnt  Irom  Report) 


i«.  surrlemcntary  notes 

♦Present  address;  Laboratory  for  Plasma  and  Fusion  Energy  Studies, 
University  of  Maryland,  College  Park,  MD  20742 
This  work  was  supported  by  the  Office  of  Naval  Research. 


1*.  key  WORDS  fConl/nu,  ofi  rfrmrmt  II  ntfttmrr  And  lamnllly  iy  block  numboi) 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


3.  RECIPIENT'S  catalog  NUMBER 


5.  TYPE  OF  REPORT  A  PERIOD  COVERED 

Interim  report  on  a  continuing 
NRL  problem. 

S.  PERFORMING  ORG.  REPORT^ DUMBER 


•  ■  contract  or  grant  NUMBERr*) 


10.  program  element,  project.  TASK 
AREA  k  WORK  UNIT  NUMBERS 

61153N;  RROll-09-41; 
47-0911-0-3 


12.  report  date 

August  24,  1983 

13.  number  of  pages 

31 _ 


IS.  security  Class.  (ol  l/U*  rmpott) 

UNCLASSIFIED 


ISa.  OECL  ASSI  Fic  ATIOH/  OOWNCRAOING 

schedule 


Level  shifts 
H 


Plasma 

Polarization 


20.  abstract  fConlinw*  on  tide  11  nmceccsfy  md  tdonttty  by  block  nvmber) 


Level  shifts  of  hydrogen  and  ionized  helium  due  to  plasma  interaction  have  been 
calculated  for  the  electron  density  10^^  cm'^  using  the  spherically  symmetrical 
distribution  of  charged  particles  around  the  radiating  atom.  The  plasma  electron 
density  has  been  obtained  from  the  solution  of  the  Schrb'dinger  equation  and  the 
ion  density  from  the  Boltzmann  distribution  function.  In  both  cases,  the  effective 


DO  1473  EDITION  OF  I  NOV  SB  IS  OBSOLETE 

S/N  0102*OI4-S601 


(Continues) 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  fWhAP  Oa(a  IntArAdJ 


1 


20  ABSTRACT  (Continued) 

potential  for  plasma  particles  has  been  determined  from  the  nuclear  charge 
of  the  radiator  and  from  the  screening  by  the  bound  electron.  Relative 
level  shifts  (i.e.,  shifts  with  respect  to  the  ground  level)  are  presented  as 
a  function  of  temperature.  The  effect  of  replacing  the  screened  potential 
by  the  true  interatomic  potential  on  the  ion  distribution  and  level  shifts 
is  discussed. 
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PLASMA  EFFECTS  ON  LEVEL  SHIFTS  OF  H  AND  He+ 
I.  Introduction 


The  interaction  of  plasma  microfields  with  highly  charged  ions  is 

attracting  considerable  attention  in  dense  plasma  physics  research, 

particularly  in  the  diagnostic  area  of  laser  driven  pellet  implosions  where 

it  is  anticipated  that  electron  densities  will  exceed  10^^  cm”^.  It  is 

fairly  well  established  that  as  the  electron  density  increases  to  values  in 

16  “3 

excess  of  say  above  10  cm  ,  a  host  of  effects  onset  that  cause  atomic 

energy  level  splitting,  shifting  and  broadening,  spectral  line  profile 
modifications,  level  merging,  and  ionization  lowering  as  well  as  a  variety 
of  other  many-body  phenomena.  All  these  effects  are,  of  course,  related 
and  probably  derive  from  a  single  unified  theory.  However,  in  the  absence 
of  such  a  theory,  each  process  is  generally  treated  individually  and 
compared  with  experimental  results  and  other  theoretical  calculations  in 
the  hope  of  gaining  further  insight  into  the  dense  plasma  physics  influence 
on  fundamental  atomic  structure  and  processes.  In  addition  to  the  basic 
and  fundamental  nature  of  this  problem,  there  is  an  immediate  and  applied 
aspect  also.  The  spectroscopy  of  dense  hot  plasmas  still  remains  one  of 
the  most  viable,  and  oftentimes  the  only,  means  of  diagnosing  conditions 
within  the  plasmas'  interior.  Therefore,  it  is  essential  that  we  explore 
the  possibilities  of  extending  conventional  theoretical  spectroscopic 
techniques  into  regions  of  higher  density. 

In  this  paper  we  will  focus  attention  on  level  shifts.  While  the 
broadening  of  spectral  lines  is  thought  to  be  better  understood  concept¬ 
ually,  the  line  shift  still  remains  a  rather  controversial  subject  both 
theoretically  as  well  as  experimentally.  We  will  assume  that  in  dense 
plasmas,  atomic  level  shifts  are  the  results  of  the  following  effects:  (a) 
Quadratic  Stark  effect,  (b)  ion  quadrupole  interaction,  (c)  electron  colli¬ 
sions,  and  (d)  plasma  polarization  in  the  vicinity  of  the  radiating  atom. 

Oftentimes  processes  (a)-(d)  collectively  are  referred  to  as  the  "plasma 

24  -3 

polarization  shift”  (PPS).  In  very  dense  plasmas  (N^  >10  cm  ),  the 

Individual  effects  cannot,  strictly  speaking,  be  separated  from  each  other, 

18  “"3 

as  is  done  at  moderate  and  low  densities  (N^  <  iu  era  ),  and  studied 
Manuscript  approved  June  28,  1983. 
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separately.  We  will  be  more  restrictive  and  assume  the  PPS  to  refer  only 
to  process  (d),  i.e.,  the  contribution  to  level  shifts  arising  from 
nonuniform  charge  distribution  in  the  neighborhood  of  the  radiator. 

The  status  of  the  PPS  prior  to  1980  has  been  reviewed  by  Griem^  and 

9 

Volonte".  A  number  of  estimates  of  level  shifts  due  to  plasma  polarization 
are  based  on  the  hydrogenic  effective  charge  approximation  in  which  the 
effects  of  electric  charges  outside  the  bound  electron  orbit  are  ignored. 
The  PPS  is  then  derived  under  various  assumptions  according  to  the  form  of 
the  free  electron  charge  distribution  inside  the  orbit.  More  recently, 
there  is  a  paper  by  Cauble  that  includes  the  electric  charge  both  inside 
and  outside  the  bound  electron  orbit.  Although  this  calculation  is  not 
self-consistent  in  the  sense  that  it  is  a  static  model  and  does  not  account 
for  charge  redistribution,  it  does  suggest  that  the  contribution  to  the 
level  shift  is  dominated  by  the  charges  outside  the  orbit  causing  a  red 
rather  than  blue  shift  of  the  L3nnan  lines  of  ionized  helium.  Also,  it  is 
worth  noting  that,  except  under  special  conditions,  the  effect  of  external 
charges  may  be  quite  different  when  the  atom  is  in  either  the  upper  and 
lower  level,  respectively,  and  may  contribute  to  the  observed  line  shift. 
A  completely  different  model  of  the  PPS  based  on  the  concept  of  the  quantum 
energy  pressure  has  been  introduced  by  Henryk,  but  even  this  model  does  not 
represent  a  fully  quantum-mechanical  solution  of  the  problem.  However, 
despite  a  number  of  simplifying  assumptions  in  ex? sting  calculations,  some 
of  the  theoretical  estimates  of  the  PPS  have  met  with  a  certain  degree  of 
success  in  the  case  of  ionized  helium  lines. 

In  this  investigation  we  will  focus  on  the  sensitivity  of  the  PPS  to 
various  forms  of  the  electron  and  ion  distributions,  particularly  on  low 
lying  levels  (n  =  1  to  4)  of  hydrogen  and  ionized  helium.  In  our  approach, 
the  energy  matrix  element  is  evaluated  for  the  interaction  of  the  single 
bound  electron  with  the  time  averaged  potential  due  to  the  plasma  electrons 
and  ions.  For  simplicity,  it  is  assumed  that  the  distribution  of  electric 
charges  around  the  radiating  atom  is  spherically  symmetric.  Similar 

1  Q  Q 

assumptions  were  invoked  by  Griem  ,  Berg  et  al.  ,  Greig  et  al.  ,  Burgess 
and  Peacock^*^,  Volonte^,  and  Pittman  et  al.^^  Both  space  and  time 
averaging  of  the  potential  is  a  questionable  procedure  for  several 
reasons.  While  the  deBroglie  wavelength  of  thermal  electrons  may  be 
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comparable  to  the  dimensions  of  the  radiator,  providing  some  justification 
for  the  continuous  nature  of  the  distribution  of  electric  charge,  the 
discrete  nature  of  the  ion  charges  is  completely  ignored  in  our 
calculation.  At  small  distances,  the  radiator  and  perturbing  ion  may 
instantaneously  form  a  quasimolecular  system  giving  rise  to  close  satellite 
lines  instead  of  a  shift  of  the  whole  line.  Also,  if  the  duration  of  the 
interaction  is  short  compared  to  the  period  of  the  orbital  motion,  then  the 
time  averaging  of  the  perturbing  potential  is  clearly  an  inadequate 
procedure.  In  spice  of  these  obstacles,  it  should  be  possible  within  the 
present  calculational  framework,  to  obtain  reasonable  estimates  of  the 
??S.  It  may  also  provide  guidance  to  more  detailed  and  mature 
investigations  in  the  future. 

18  "“3 

In  low  density  plasmas  (N^  <10  cm  )  the  Debye  shielding  distance 
is  much  larger  than  the  dimensions  of  the  radiating  atom  and  the  level 
shifts  are  produced  mainly  by  the  charge  distribution  in  the  vicinity  of 
the  radiator.  The  method  utilized  here  treats  the  motion  of  the  plasma 
electrons  quantum  mechanically  while  the  ions  are  described  by  classical 
Boltzmann  statistics.  Also,  the  potential  for  both  the  free  electrons  and 
ions  is  derived  from  the  nuclear  charge  and  from  the  screening  by  the  bound 
electron.  The  effects  of  the  ion  distribution  based  on  interatomic  poten¬ 
tials  are  included  in  several  examples.  Finally,  Che  mutual  Interaction  of 
plasma  electrons  as  well  as  that  of  the  ions  is  neglected  and  consequently 
near  the  radiator  no  account  is  taken  of  collisional  redistribution  of 
electron  velocities. 

II.  General  Considerations 

In  this  calculation  of  atomic  level  shifts  in  one-electron  systems  due 
to  the  plasma  polarization  effect  we  adopt  the  simplifying  assumption  that 
the  actual  distribution  of  perturbing  charges  around  the  radiating  atom  can 
be  replaced  by  a  continuous,  spherically  S5mimetric  time  independent  charge 
distribution.  Even  with  this  assumption,  the  problem  of  calculating  atomic 
level  shifts  may  be  approached  from  two  different  viewpoints:  either  the 
charge  distribution  around  the  radiator  corresponds  to  an  average  state  of 
excitation  of  the  radiating  atom  (case  A),  or  it  corresponds  to  a  partic¬ 
ular  level  n£  of  the  radiator  (case  B).  Case  A  is  equivalent  to  the 
assumption  that  plasma  particles  do  not  have  time  to  adjust  to  Individual 
atomic  orbitals. 
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We  will  assume  that  only  singly  charged  ions  are  present  in  the  plasma 
which  may  be  partially  ionized,  and  that  the  effect  of  neutral  atoms 
surrounding  the  radiator  may  be  neglected. 

The  total  energy  of  the  system  consisting  of  the  nji-bound 
electron  in  the  field  of  the  nuclear  charge  Z  and  the  surrounding  plasma  in 
a  given  volume  may  be  written  as  a  sum 


E  =  E,  +  E-  +  E.  +  E.  , 
njl  1  2  3  4 


(1) 


where  E^  is  the  expectation  value  of  the  operator  -  y  V"  -  Z  r  (in 
atomic  units,  which  are  used  throughout  the  paper,  unless  indicated 
otherwise)  evaluated  with  the  wave  function  of  the  bound  electron  in 
the  n£  orbit,  E2  represents  the  electrostatic  interaction  between  the  bound 
electron  and  the  plasma,  E^  is  the  interaction  energy  of  the  nucleus  and 
the  plasma,  and  is  the  sum  of  all  contributions  not  explicitly  included 
in  E^,  E2,  or  E^.  Effects  of  exchange  between  the  bound  and  free  electrons 
will  be  omitted  in  the  present  calculations.  In  the  first  approximation, 
the  energy  of  an  emitted  or  absorbed  photon  is  equal  to 


("1 


+  E2  +  E3  + 


(2) 


where  the  quantities  in  parentheses  (  and  (  ) to  be  evaluated 
with  respect  to  the  nl  and  n'Z'  atomic  level,  respectively.  In  most 
experimental  plasma  conditions  in  which  neutral  hydrogen  and  singly  ionized 
helium  are  observed,  the  perturbation  of  atomic  wave  functions  is  very 
small  and  we  can  assume  that  (E.)  is  equal  to  the  eigenvalue  e  of  the 

It  L  ^  ^  ^ 

unperturbed  Schrodinger  equation  for  the  hydrogen-like  atom. 

The  frequency  shift  Av  of  a  spectral  line  due  to  the  plasma  inter¬ 
action  is  then  equal  to 


Av 


h"^(E  -E 
njl 


n  I' 


+  e  )  . 

n'l' 
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Because  only  energy  differences  are  directly  observable,  it  is  convenient 

to  express  all  energies  E  with  respect  to  the  energy  of  the  ground  level 

n£ 

Is.  As  a  relative  shift  of  the  aZ  level  we  define  a  quantity 


(3) 


so  that 

Av  =  h'^  ae;^  . 

if  njl>ls  is  an  optical  transition. 

III.  Method  of  Calculation 


Define 

W(r)  *  Zr  ^  +  4Tir  p(r')  r'^dr  +  p  (r  )  r  dr'  .  (4) 

p(r')  =  where  and  p^  are  the  spherically  averaged 

charge  densities  of  the  plasma  electrons  and  ions,  respectively.  The 
potential  energy  of  the  bound  electron  in  the  field  of  plasma  particles 
is  given  by 

V.  =>  -W  +  Zr~^  .  (5) 

D 

We  assume  that  the  plasma  electrons  have  a  Maxwellian  velocity 
distribution  at  large  distances  from  the  nucleus  and  that  in  the  vicinity 
of  the  radiating  atom,  in  a  low-density  plasma,  p^  and  p^  are  determined 
only  by  the  electrostatic  potential  of  the  nucleus  and  the  bound 
electron.  The  mutual  interaction  of  plasma  electrons  and  ions  will  be 
Ignored.  Consequently,  the  potential  used  for  the  calculation  of  p^  and  p^^ 
behaves  like  (Z-l)/r  at  large  r,  while  in  a  real  plasma  this  interaction 
leads  to  the  formation  of  a  charged  cloud  and  to  a  faster  decrease  of 
potential.  In  low-density  helium  plasmas  with  Z=2  the  dimensions  of  this 
cloud  are  several  orders  of  magnitude  larger  than  the  size  of  low  lying 
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atomic  orbitals  and  the  effect  of  distant  charges  on  relative  atomic  energy 
levels  will  be  ignored. 

In  our  procedure,  is  defined  by 

Pg(r)  =  f(k)  (kr")  ^  (2£+l)  dk,  (6) 

where  f(k)  is  the  velocity  distribution  of  free  electrons  normalized  to 

unity,  Ng  is  the  average  number  density  of  plasma  electrons,  and  the 

functions  F,  are  solutions  of  the  equation 
R  % 

_  ILW  -2V^  +  k^]  F  (r)  -  0  (7) 

dr'^  r'^ 

with  the  asymptotic  form 

FkJl(r)  ^  sin(kr+6^^)  .  (8) 

For  scattering  on  a  neutral  target,  is  independent  of  r.  From  (6),  (7) 

and  (3)  it  follows  that  p^>  at  large  values  of  r.  in  Eq.  (7)  is  the 
potential  energy  of  the  plasma  electron  in  the  field  of  the  nucleus  and  the 
hound  electron.  If  p^  is  the  spherically  averaged  charge  density  of  the 
bound  electron,  is  given  by 

Vg(r)  =  -W(r)  (9) 

with  p  in  expression  (4)  for  W(r)  replaced  by  p^. 

CASE  A. 

In  this  case,  and  are  independent  of  nji  and  AE'  is  given  by 

11  Xr 

P  (r)  and  P,  (r)  are  unperturbed  radial  functions  of  the  bound  electron 
n  x.  -i-S  2 

normalized  to  satisfy  the  condition  J  P^^  dr  =  1,  and  V^(r),  given  by  Eq. 
(5),  is  the  potential  energy  of  the  bound  electron  at  a  distance  r  from  the 
nucleus  in  the  field  created  by  plasma  electrons  and  ions. 
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p,  ::iay  be  written  as 

D 

-  (4rr")~^  Z  b  (r),  (IL) 

ni  ni,  nl 

where  the  occupation  numbers  satisfy  the  condition  Zb  =1.  In  a  low- 

nZ  nJl 

density,  low-temperature  plasma,  the  occupation  of  excited  states  is  very 
small  and  we  will  assume  that  only  the  ground  level  Is  is  occupied. 

With  the  assumption  that  only  singly  ionized  ions  are  present  in  the 
plasma,  we  adopt  the  following  form  for  the  positive  ion  charge 
distribution: 


p.  =  exp  (-V./kT),  (12) 

where  V^(r)  is  the  interatomic  potential  for  the  radiating  atom  and  a 
plasma  ion.  In  our  procedure 


V^(r)  =  -V^(r). 


(13) 


CASE  B. 


All  three  quantities  Ep,  E^  and  E^  may  depend  on  nz  as  well  as  the 
potential  energy  Vjj(r).  We  restrict  our  investigation  to  the  calculation 
of  the  quantity  (l'2^nZ  ”  '^2^1s  case  A,  and  adopt  the  assumption 

that  the  sum  E^+E^  does  not  depend  on  nz.  Then 

„  =  /  [(V,)  -  (V,),  pj  ]  dr  (14) 

nz  •'  b'nZ  nZ  b  Is  Is'*  ^  ^ 

with  defined  by  (5).  Furthermore, 

-  -  (4^r^)"^  (r)  ,  (15) 

and  p^,  are  defined  by  (12)  and  (13).  In  case  B,  both  and  depend 
on  the  state  of  the  radiator. 
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Eq.  (13)  implies  that  the  bound  state  wave  function  P  remains 

nz 

essentially  unperturbed  by  plasma  ions.  Other  forms  of  interatomic 
potential  will  be  discussed  in  the  Appendix. 

The  charge  profile  in  (6)  was  obtained  numerically  from  7  different 
values  of  k.  for  which  Eq .  (7)  has  been  solved.  Number  of  angular  momenta 
included  in  the  summation  over  £  in  (6)  varied  from  10  to  50  for  hydrogen 
and  from  20  to  60  for  ionized  helium,  depending  on  tc. 

IV.  Results  for  Hydrogen  and  Ionized  Helium. 


Calculations  of  relative  level  shifts  have  been  made  for  hydrogen  for 

temperatures  from  0.5  to  2.0  eV  and  for  singly  ionized  helium  from  2.5  to 

4.0  eV.  Since  the  relative  shifts  in  low-density  plasmas  are  proportional 

to  the  average  electron  density  N^  (see  Eqs.  ,4),  (5),  (10),  (12)),  our 

17  — q 

results  are  normalized  to  =  10  cm  .  Relative  level  shifts  AE"  for 

®  n£ 

case  A  are  shown  in  Figs.  1  and  3  for  H  and  He"*”,  respectively,  and  for  case 
B  are  displayed  on  Figs.  2  and  4. 

Plasma  interactions  remove  the  degeneracy  of  £-sublevels.  In  our 
approximation  all  effects  of  level  broadening  are  ignored  and  the  "plittlng 
is  well  defined. 

In  case  A,  relative  level  shifts  for  both  hydrogen  and  ionized  helium 

are  negative  and  lead  to  the  red  shift  of  spectral  lines  of  the  Lyman 

series.  The  i-splitting  of  relative  level  shifts  due  to  plasma 

interactions  is  generally  larger  for  He'*'  than  for  H,  and  the  values  of 

AE' ^  increase  with  £. 
n£ 

In  case  B,  the  behavior  of  relative  shifts  for  hydrogen  is  drastically 
changed.  Almost  all  values  of  are  positive,  the  ^.-splitting  is  so 

large  that  some  curves  corresponding  to  different  principal  quantum  numbers 
cross  each  other,  and  relative  shifts  AE'„  increase  with  decreasing  £.  For 

nit 

ionized  helium,  case  B  is  not  much  different  from  case  A  due  to  the 

dominant  attractive  force  of  the  nuclear  charge.  All  aE'  values  are 

n£ 

slightly  shifted  upwards  but  remain  negative  as  in  case  A. 

When  comparing  the  present  results  with  the  observed  line  shifts,  one 
has  to  bear  in  mind  that  (a)  using  differences  of  lev^l  shifts  to  calculate 
line  shifts  represents  only  an  approxi.^ation,  (b)  other  effects  mentioned 
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in  the  introduction  contribute  to  the  line  shifts  and  may  be  larger  or 

comparable  to  the  values  of  AE'  ,  and  (c)  our  relative  shifts  aE'  in  case 

n£’  nz 

B  (Figs.  2  and  4)  do  not  include  contributions  from  E^  and  E^  in  Eq.  (2). 

For  He"*”,  our  results  predict  a  very  small  red  shift  of  Lyman  series 

lines.  The  shifts  are  of  the  same  order  of  magnitude  as  those  obtained 

from  the  Debye-Huckel  potential.  The  experimental  evidence,  however, 

points  to  much  larger  blue  shifts  ’  .  For  the  3-2  transition  of  He  at 

1640  4  our  calculations  predict  a  negligible  red  shift  in  agreement  with 

experiment^,  but  the  theoretical  shift  for  the  4-2  transition  at  1215  A 

is  again  too  small,  except  possibly  for  the  temperature  around  4  eV. 

Pittman  et  al.^^’^^  measured  a  red  shift  of  0.17  A  for  the  transition  4-3 

at  4686  A  for  T  *  4  eV  in  He"*".  However,  the  ion  contribution  itself 

(quadrupole  and  quadratic)  has  been  estimated  to  be  -0.24  A,  so  that  our 

theoretical  result  +0.11A  for  the  4f-3d  component  (Fig.  3)  for  the  PPS 

cannot  explain  the  observed  red  shift.  Wiese  et  al.^^  also  measured  red 

shift  of  Balmer  series  lines  of  hydrogen.  Our  theoretical  results  predict 

blue  shifts  for  H  and  either  blue  or  red  shifts  for  Hg,  depending  on 

a 

temperature  and  line  component,  but  the  calculated  values  are  much  smaller 

than  the  observed  ones.  Measurements  of  the  L  shifts  of  hydrogen  have 

13  ^ 

been  reported  by  Fussman  .  The  blue  shift,  if  any,  was  found  to  be  less 
than  0.015A  at  T  ■  leV  which  agrees  with  our  theoretical  estimates. 

We  arrive  at  the  conclusion  that  the  PPS  for  hydrogen  and  ionized 
helium  as  defined  in  the  present  investigation  is  too  small  to  be 
responsible  for  the  majority  of  measured  line  shifts  in  contrast  to  the 
results  of  theories  based  on  the  hydrogenic  effective  charge  approxi¬ 
mation,  and  that  other  mechanisms  have  to  be  taken  into  account  in  order  to 

obtain  an  agreement  with  observations.  A  similar  conclusion  has  been 
3  + 

reached  by  Cauble  who  calculated  the  PPS  for  He  by  a  method  similar  to 
that  of  Volont^^,  but  included  screening  effects  of  plasma  electrons 
outside  the  orbit  of  the  bound  electron. 

In  the  present  investigation,  the  potential  energy  for  the  ion 
distribution  in  (12)  was  derived  from  the  screening  by  the  bound  electron, 
and  it  was  assumed  that  the  bound  electron  wave  function  remained 
unperturbed  during  collisions  with  ions.  Consequently  there  was  no 
difference  between  level  shifts  obtained  for  helium  ions  in  a  pure  helium 
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plasma,  and  for  helium  ions  in  a  hydrogen  plasma  containing  helium  ions 
only  as  impurities.  In  the  Appendix  we  will  investigate  the  effects  of 
ions  on  level  shifts  by  substituting  actual  interatomic  potentials  for 
in  Eq.  (12). 
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V.  Appendix 


During  close  collisions  with  ions,  the  bound  electron  wave  functions 
may  be  appreciably  distorted  by  the  strong  electric  field  of  the  colliding 
ion.  In  some  cases  a  quasimolecule  may  be  temporarily  formed.  The  short- 
range  effects  may  influence  the  distribution  of  ions  around  the  radiating 
atom  and  can  be  taken  into  account,  in  the  first  approximation,  by 
replacing  in  Eq.  (12)  by  true  interatomic  potentials. 


Generally,  more  than  one  potential  curve  corresponds  to  a  particular  ni 
state  of  a  hydrogen  atom  or  a  hydrogenic  ion  and  a  plasma  ion  in  its 
ground  state.  Each  of  these  potential  curves  is  associated  with  a 
different  state  of  the  united  atom  as  the  nuclear  separation  approaches 
zero.  E.g.,  there  are  two  potential  curves  correspponding  to  a  hydrogen 
atom  in  the  ground  state,  and  a  proton,  associated  with  the  and  2pa^ 

states  of  the  molecule.  Instead  of  using  Eq.  (12),  in  this  case  we 
calculate  the  ion  charge  distribution  from  the  relation 


Pi  “  Ng(5  5  gj  exp  [-V^(A^)/kT]  ,  (16) 

where  ^^(Aj)  :*s  the  interatomic  potential  energy  corresponding  to  the 

molecular  state  A^  with  the  statistical  weight  gj,  and  the  summation  is 
carried  out  over  all  molecular  states  leading  to  the  same  state  of 
separated  atoms. 

Potential  curves  Vi(Aj)  for  helium  ions  in  a  pure  helium  plasma 
correspond  to  the  process  He2  He  (njl)  +  He  (Is).  However,  sufficient 
data  on  interatomic  potentials  for  excited  states  of  He2  are  not 
available,  and  therefore  we  restrict  our  calculations  of  He'*'  level  shifts 
to  the  case  of  impruity  helium  ions  in  a  hydrogen  plasma.  In  this  case  we 

44  4  4 

consider  potential  curves  for  the  process  HeH  -*•  He  (nZ)  +  H  ,  and  use 

potentials  calculated  by  Bates  and  Carson^^  for  ni  =  Is,  2s,  and  2p. 

For  relative  level  shifts  in  hydrogen,  we  consider  the  process  H^ 

->•  H(n)l)  +  H  .  The  interatomic  potentials  are  taken  from  Bates  et  al.^^  for 

ni  *  la,  2s,  2p,  and  3s.  The  published  potential  curves  were  extrapolated 
to  larger  values  of  nuclear  separation,  where  needed. 


The  interatomic  potentials  V^(Aj)  are  compared  with  the  screened 
potentials  used  in  par.  4  in  Figs.  5  -  10. 

The  ion  distribution  calculated  from  (16)  using  interatomic  potentials 
has  a  profound  effect  on  the  relative  level  shifts.  As  in  the  text,  we 
consider  two  cases:  in  case  A,  the  relative  level  shifts  AE',  for  all  nA 
levels  are  calculated  from  the  same  V^{r).  The  distribution  of  plasma 
electrons,  given  by  (6),  is  obtained  from  the  screening  by  the  Is  bound 
electron,  and  the  distribution  of  ions,  given  by  (16),  is  obtained  from  the 
interatomic  potentials  corresponding  to  H(ls)  +  H  or  He  (Is)  +  H  , 
respectively.  In  case  B,  the  potential  energies  V^(r)  corresponding  to 
different  ni  levels  are  different.  For  a  given  nl  level,  the  distribution 
of  plasma  electrons  is  determined  from  the  screening  by  the  n£,  bound 
electron,  and  the  ion  distribution  by  interatomic  potentials  corresponding 
to  H(n)l)  +  h"*”  or  He'^(njl)  +  h"*".  Results  for  AE'  calculated  from  (10)  or 
(14),  respectively,  are  shown  in  Figs.  11-14. 


The  attractive  character  of  the  potential  curve  for  the  Isa  state  of 
j.  S 

H2  at  small  r  (Fig.  5)  leads  to  the  accumulation  of  protons  in  the  vicinity 

of  the  hydrogen  atom  in  the  ground  state  and  to  the  lowering  of  the  Is 


level. 


Therefore  all  relative  shifts  AE 


ni 


rapidly  increase  with 


decreasing  temperature  as  shown  on  Figs.  11  and  12.  The  large  difference 

of  aEI  and  aEI  in  hydrogen  (case  B,  Fig.  12)  is  the  result  of  the  same 
is  Zp 

effect;  potential  curves  of  H(2s)4H'’’  are  repulsive,  while  the  curves 

states  of  the  H(2p)+H'''  system  exhibit 


representing  the  2pir  and  3da 

u  g 

minima.  Consequently  the  upward  shift  of  the  2p  level  is  reduced  by  the 
accumulation  of  positive  charges,  and  AE'^  <  AE'^. 

The  potential  curve  corresponding  to  He'*’(ls)-ttl'*’  has  no  minimum  and  the 
values  of  relative  shifts  for  He'*’  (Figs.  13,14)  do  not  increase  towards 
lower  temperature  in  contrast  to  hydrogen.  This  interatomic  potential  is 
very  similar  to  the  potential  obtained  from  the  simple  screening  by  Is 
electron  (Fig.  8)  and  therefore  there  is  not  much  difference  between  AE'^ 
on  Fig.  3  and  13.  The  large  differences  of  AE'^  and  AE'p  in  case  B  (Fig. 
14)  again  arises  from  the  difference  of  potential  curves  associated  with 
the  He^(2s)+ft^  and  He^(2p)+H'^  systems. 
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I”  spite  of  these  changes  in  relative  level  shifts,  the  general 
conclusion  about  the  magnitude  of  the  shifts  obtained  within  our 
approximation,  which  has  been  discussed  earlier,  remains  valid.  The 
calculated  relative  shifts  of  levels  for  which  the  interatomic  potentials 
are  known,  appear  to  be  small  and  additional  mechanisms  must  be 
considered  for  the  explanation  of  measured  line  shifts. 
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kT  (eV) 

Fig.  1.  Relative  level  shifts  AE'  for  hydrogen.  Hydrogen  plasma,  N  = 

17  -3 

10  cm  .  Case  A  (All  relative  shifts  are  calculated  from  the  same  perturbing 
potential.  The  distribution  of  plasma  electrons  and  ions  around  the  radiating 
atom  is  determined  from  screening  by  the  Is  bound  electron).  The  curves  for 
4p  and  4d  lie  between  those  for  4s  and  4f.  They  are  not  shown  for  clarity. 
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0-  kT(eV) 

,  1  for  hydrogen.  Hydrogen  plasma,  - 

?lg.  2.  Relative  level  shift  j  calculated  free  a 

Case  B  (Each  quantity  CEg)^  derived  fre. 

potential  corresponding  to  the  electron  a 
screening  by  the  nl  bound  electron). 


RELATIVE  LEVEL  SHIFT  (10"® RYD) 


0 


Fig.  3.  Relative  level  shifts  for  ionized  helium, 

plasma.  -  lO^^cm”^.  Case  A. 


17 


Hydrogen  or  helium 


rV|(r) 


Fig.  5.  Effective  charge  acting  on  the  proton  in  the  system  H(ls)  +  Full 

curves;  Interatomic  potentials  from  Ref.  15.  Dashed  curve:  Potential 
derived  from  screening  by  the  Is  bound  electron. 
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Fig.  6.  Effective  charge  acting  on  the  proton  in  the  systems  H(2s)  + 

Full  curves:  from  Ref.  15.  Dashed  curve:  from  screening  by  the  2s 

electron. 
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rV,(r) 


I 


Fig.  7.  Effective  charge  acting  on  the  proton  in  the  system  H(2p)  +  H"*”.  Full 
curves;  from  Ref.  15.  Dashed  curve:  from  screening  by  the  2p 

electron. 
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Fig.  8.  Effective  charge  acting  on  the  proton  in  the  system  He""(ls)  +  H' 
Full  curve:  from  Ref.  14.  Dashed  curve:  Vj_  from  screening  by  the  Is 
electron. 


rV,(r) 


Fig.  10.  Effective  charge  acting  on  the  proton  in  the  system  He'*’{2p)  +  H'*’. 
Full  curves:  from  Ref.  14.  Dashed  curve:  from  screening  by  the  2p 
electron. 
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Fig.  IL.  Relative  level  shifts  AE''  for  hydrogen.  Hydrogen  plasma,  N_  * 

17  —3 

10  cm  .  Case  A  (Same  as  in  Fig.  1,  but  the  proton  distribution  is  derived 
from  interatomic  potentials  corresponding  to  H(l8)  +  H'*'  shown  in  Fig.  5).  T 
curve  for  3p  lies  between  those  for  3s  and  3d,  and  the  curves  for  4p  and  4d 
lie  between  48  and  4f.  They  are  not  shown  for  clarity. 


Fig.  12.  Relative  level  shifts  AE'  for  hydrogen.  Hydrogen  plasma,  N  = 

17  _3 

10  cm  .  Case  B  (Same  as  in  Fig.  2,  but  the  proton  distribution  is  derived 
from  interatomic  potentials  corresponding  to  H(2s)  +  H''",  H(2p)  +  and  H(3s) 
+  H"*",  respectively). 


kT(eV) 

Fig.  13.  Relative  level  shifts  for  ionized  helium  in  a  hydrogen  plasma, 

17-3 

Ng  •  10  cm  .  Case  A  (Same  as  in  Fig.  1,  but  the  proton  distribution  is 
derived  from  the  interatomic  potential  corresponding  to  He'*’(ls)  +  H"*”  shown  in 
Fig.  8). 
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Fig.  14.  Relative  level  shifts  AE'  for  Ionized  helium  in  a  hydrogen  plasma, 
17  —3 

Ng  “  10  cm  .  Case  B  (same  as  in  Fig.  2,  but  the  proton  distribution  is 
derived  from  interatomic  potentials  corresponding  to  He  (2s)  +  H  and  He  (2p) 
+  H’'",  respectively,  shown  in  Figs.  9  and  10). 


